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ABSTRACT 

We describe a pseudo-Newtonian potential which, to within 1% error at all angular momenta, 
reproduces the precession due to general relativity of particles whose specific orbital energy is small 
compared to in the Schwarzschild metric. For bound orbits the constraint of low energy is equivalent 
to requiring the apoapsis of a particle to be large compared to the Schwarzschild radius. Such low 
energy orbits are ubiquitous close to supermassive black holes in galactic nuclei, but the potential 
is relevant in any context containing particles on low energy orbits. Like the more complex post- 
Newtonian expressions, the potential correctly reproduces the precession in the far-field, but also 
correctly reproduces the position and magnitude of the logarithmic divergence in precession for low 
angular momentum orbits. An additional advantage lies in its simplicity, both in computation and 
implementation. We also provide two simpler, but less accurate potentials, for cases where orbits 
always remain at large angular momenta, or when the extra accuracy is not needed. In all of the 
presented cases the accuracy in precession in low energy orbits exceeds that of the well known potential 
of Paczyhski & Wiita (1980), which has ^ 30% error in the precession at all angular momenta. 

Subject headings: Methods: numerical — black hole physics — galaxies: kinematics and dynamics — 
galaxies: nuclei 



1. INTRODUCTION 

Pseudo-Newtonian potentials that modify the Newto- 
nian gravitational potential have a long history of use in 
astrophysics. While general relativity is now well under- 
stood in the astrophysics community, pseudo-Newtonian 
potentials are still useful in approximating relativistic ef- 
fects in sim pler and faster Newtonian simulations. The 
potential of iPaczvnski fc Wiital (|1980( ) is often used in 
the study of accretion onto relativistic objects. In this 
regime the Paczyhski- Wiita potential often gives results 
close to those using full GR since it correctly repro- 
duces the location of the inner most stable circular orbit 
(ISCO) and the marginally bound orbit as well as being 
a good approxim ation to the bindin g energy at the ISCO 
(for a review see I Abramowicd 120091 ) . 

Here we propose a series of Newtonian potentials with 
a different aim: to correctly reproduce the precession 
produced by general relativity in the Schwarzschild met- 
ric for test particles whose apoapsis lies far from the 
hole, i.e., in the non-relativistic region. The Paczynski- 
Wiita pot ential has been u sed in this context multiple 
times (e.g. lChen et al.ir2011[ ). despite its key properties of 
closely reproducing the location and energy of the ISCO 
being unimportant in this regime. Instead, we propose 
alternative potentials that are more accurate and physi- 
cally relevant for these orbits. We have used them to sim- 
ulate galacti c dynamics around su permassive black holes 
(SMBHs) in lBode fc Weggl(|20ll . The primary concern 
in that context was to ensure that stars passing close 
to the black hole exited along the correct trajectories. 
These potentials are likely to be useful in other contexts, 
motivating the brief presentation here. Throughout we 
use geometrized units where G ^ c = 1 . 

2. SUMMARY OF PROPOSED POTENTIALS 
^ [wegg@tapir.caltech.edu| 
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Figure 1. Comparison of the precession per orbit produced by the 
proposed potentials with the GR expression for parabolic orbits as 
a function of specific angular momentum, h. The labeled potentials 
are described in Table m The GR precession is shown by a solid 
black line (in the upper panel the precession produced by potential 
C lies almost on top of the GR expression). In the lower panel 
we plot the fractional error relative to the relativistic precession, 
defined to be (A(^ - A(/)GR)/(A(^Ga - 27r). 



We present three new pseudo-Newtonian potentials in 
this paper. All of these potentials, and the potential of 
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Table 1 

Coefficients for the potentials U{r) = —aAI/r — (1 — a)M /(r — Rx) — MRy/r^ {equation[TJ described in this work. 

Potential a R^/IM Ry/M Precession error: r 3> M Diverges at = AM Maximum precession error 
A i ~ 3 No 100%=^ 
B 5/3 4/3 Yes < 30% 
C ^ (2 + Ve) (4^6 - 9) ^ - 3) Yes < 1% 
Paczyhski-Wiita 2 33% Yes 33% 



''Does not diverge at h = 4M 



IPaczvnski fc Wiital (fT980t) . can be written in the form 
aM (1 - a)M MRy 



U{t) 



(1) 



where the values of the coefficients a, Rx and Ry for 
the potentials are summarized in Table [H We choose 
potentials of this form since the presence of the 
term allows the far field behavior to reproduced, while 
the l/(r — Rx) term allows reproduction of the divergent 
behavior as the specific angular momentum approaches 
4M. The resultant precession per orbit is compared to 
the GR value in Figure [T] In what follows we justify 
these choices. 

3. APPROACH TO CALCULATING PROPOSED 
POTENTIALS 

3.1. Precession Due to General Relativity 

In general relativity the change in azimuthal angle of 
a test particle between two consecutive apoapsides on 
a geodesic in the Schwarzschild metric is given by (e.g. 
Equation 25.42 of[AIisner et al. 1973.) 



A0GR 



2M 



-1/2 



dr 



(2) 

where E is the specific energy of the particle without rest 
mass energy (i.e. E = —po //i — 1 where fj. is the particles 
mass), h is the specific angular momentum (i.e. h = 
p^/ and r± are the radii of periapsis (— ) and apoapsis 
(+) given by the two largest roots of the equation 
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= 0. 



(3) 



For our 'nearly parabolic orbits' {E <C 1) the precession 
due to relativity is a function only of the angular mo- 
mentum. Unless otherwise noted in this paper we work 
in the limit that E = 0. 

3.2. Precession due to Newtonian Central Potential 

By comparison, in classical mechanics the change in 
azimuthal angle for a test particle of any energy between 
two consecutive apoapsides in a central potential, C/(r), 
is given by (e.g. iLandau fc Lifshit3ll969D 
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where in this case r± are given by 



2{E-U{r±)) 



= 0. 



r±' 



(5) 



3.3. Requirements of Proposed Pseudo- Newtonian 
Potentials 

In principle it is possible to define a pseudo-Newtonian 
potential, U{r), such that the precession angles given by 
equations [2] and m are equal in the limit <C 1, i.e. 



A(^gr(^) = A(/)Ncwt.(/l) 



(6) 



This potential would have the property desired: on re- 
turning to large radii, test particles would have precessed 
through the correct angle, and be traveling along the cor- 
rect path with only a time error. However, we also desire 
a simple potential for efficient calculation and so instead 
we seek potentials that minimize the precession error, i50, 
defined through 



<50(/l) = A(/)Gr(/i) - A0Ncwt.(M 



(7) 



We propose three potentials that, in order of complex- 
ity, minimize 5<j): A) in the far field (large h) ; B) in the 
far field and whose precession diverges logarithmically in 
the same location as GR {h — ?> 4M); and C) in the far 
field, and whose precession diverges logarithmically as 
h AM with the correct magnitude. 

4. PROPOSED PSEUDO-NEWTONIAN 
POTENTIALS 

4.1. Potential A: Matching The Far Field Precession 

In this section we consider the behavior of orbits with 
h » 4A/, but we do not require E ~ 0, only that E <^ 1. 
In this case, inspection of Equation [3] shows r± ^ M 
and the entire orbit lies in the far field. 

The change in angle in the far field in GR can be cal- 
culated from e quation [2] and is well known to be (e.g. 
IWeinbersl[T97l) 



A0gr(/i) = 27r 



/l2 



for h > UI . 



(8) 



Note that all that is required is a sufficiently distant pe- 
riapse. It is not required that the orbit have E = 0. 
In the far field we require potentials have the form 
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M MR,, 



2 +0(r-^). (9) 
Neglecting the terms 0{r^'^) and higher, after some al- 
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gebra, the precession calculated from equation 2] is 

2r,,a 



A0Ncwt.(^) = 2 



(r — ?'_)(r+ — r) 



(10) 



where r,, = h'^/2M, a = ~M /2E, and 



r± = a 



l±\ll--{Ry-rh) 



(11) 



Contour integration gives the integral 



A</>Newt.(/i) =27r( 1 " 



27r- 



2TTMRy 



-1/2 



+ Oia/rh). 



(12) 



Provided that h ^ M and i? ^ 1 , then the final term can 
be dropped and matching the far field precession given 
by equation [5] requires Ry = 3M. 

When concerned with the far field precession we there- 
fore propose the pseudo-Newtonian potential 



C/(r) 



1 



3M 



(13) 



which requires only modest additional computation over 
the Newtonian expression. Note that this pseudo- 
Newtonian potential gives the correct precession for all 
orbits (parabolic, eccentric, or hyperbolic), provided that 
the periapse lies in the far field. The resultant precession, 
labeled as potential A, is plotted in Figure [TJ 

Since this potential does not reproduce the diver- 
gence as /i — ^ AM this potential performs worse than 
the Paczynski-Wiita potential, even for parabolic orbits, 
when h < 4.8M which corresponds to periapsis separa- 
tion of r < 8.5M. For reference h = 4M corresponds to 
periapsis separation r = 5M for parabolic orbits in this 
potential. 

4.2. Potential B: Logarithmic Divergence As h AM 

In this section we construct a potential that reproduces 
the logarithmic divergence of the general relativistic pre- 
cession as /i — 4M. 

First consider a potential of the form 

t/(.) = --iL-. (14) 
The precession in this potential for E' = is given by 

^2 1-1/2 
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dr 
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The roots of the quadratic form in the integrand can be 
written as 
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SO that 

A0Ncwt.(^) = 2 



rhjr - R^) 
{r-rp){r-rq) 
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dr 
r 



(16) 



(17) 



As ru 



ARx- then r„ 



and the integral diverges 



logarithmically, rt ARx corresponds to h ^ y/8MRx 
and the leading order behavior of the integral is 



lim A(j>-!^cwt.{h) 



log h - ^/8MR, 



(18) 



A similar calculation using the GR expression gives a 
logarithmic divergence as /i — > AM and the correspond- 
ing expression is 



lim A(f>GR{h) 



-V2log{h~AM) 



(19) 



For A(/)Newt. to diverge at h = AM, we must have R^ 



2M. T his is exactly the potential of iPaczvhski fc Wiital 
(|1980( ). which therefore diverges at the correct angular 
momentum. H owever comparing equatio ns 1181 and 1191 
the potential of IPaczvhski &: Wiital ()1980[) has the incor- 
rect magnitude (by a factor of V2) as the angular mo- 
mentum approaches AM, and as noted previously, has 
incorrect precession in the far field. 
To correct the far field behavior consider the potential 



U{r) = - 
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The calculation of the precession proceeds in the same 
manner, but with the roots now given by 
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rh - Ry 



(21) 



Again, as the roots coincide the integral diverges log- 
arithmically. Requiring that the divergence occurs as 
h — >■ AM and that the potential has the correct far 
field limit (i.e. far from the hole the expansion is given 
by equation [9] with Ry — 3Af) gives R^ = 5M /3 and 
Ry — AM /3. Our proposed potential which has the cor- 
rect precession in the far field and which also logarithmi- 
cally diverges as ft. — >■ AM is therefore 
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(22) 



The resultant precession, labelled as potential B, is 
shown in Figure [TJ This is the potential used in 
iBode fc Weg g (2011]). For reference h = AM corresponds 
to periapsis separation r ~ 10il//3 for parabolic orbits 
in this potential. 

4.3. Potential C: Correct Rate Of Logarithmic 
Divergence As h ^ AM 

The potential proposed in equation [22l has the correct 
far field behavior, and diverges logarithmically at the 
correct angular momentum. However, the rate of that 
divergence is incorrect: The behavior of the integral as 
h AM is 



lim A(/)Ncwt.(ft) 

h^iM+ 



AM) 



(23) 



which does not match the GR expression in equation ll9l 
An additional term in the potential allows this to rec- 
tified. Using a potential of the form 



U{r)^ — 



aM (1 - a)M MRy 



Rx 



(24) 
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enables us to match the three constraints for the three 
coefficients a, Rx and Ry. The constraints on the co- 
efficients are that: 1) in the far field the precession ap- 
proaches equation [8] (i.e. equation [9] with Ry — 3A/), 

2) the integral diverges logarithmically as h ^ AM, and 

3) the rate of divergence as ft. — )• 4M is given by equa- 
tion [ini The values of a, Rx and Ry satisfying these 
constraints are 



= ( 4^/6 - 9 ) M 



(25) 



Ry ^ ^ f2V6-3)M. 



The precession produced by this "potential C" is com- 
pared to the GR expression in Figure [1] For refer- 
ence h = AM corresponds to periapsis separation of 
r = 2(\/6 — 1)M for parabolic orbits in this potential. 

This potential produces precession which agrees with 
GR to within 1% for all orbits where E <^ 1 i.e. 
whose specific orbital energy is small compared to in 
the Schwarzschild metric. For bound orbits this corre- 
sponds to requiring apoapsis be large compared to the 
Schwarzschild radius (r+ ^ M) . 

5. CONCLUSIONS 

We have proposed three pseudo-Newtonian potentials 
appropriate for 'nearly parabolic orbits' (orbital energy, 
i? <C 1) around a Schwarzschild black hole. These nearly 
parabolic orbits correspond to orbits with large apoapsis 
compared to the Schwarzschild radius of a central black 
hole, or mildly hyperbolic orbits. 

For bodies which pass close to the black hole, these po- 
tentials accurately reproduce the changes in the "New- 
tonian" parts of the trajectories far from the black hole, 
differing from the exact GR expression only by a time 
error. In the far field the time error as a fraction of the 
orbital period is of order 5P/P = 0{E) < 1 but di- 
verges as /i ^ AM like 5P/P = 0{E^^^\og{h - AM)). 
Therefore, for i? <C 1, the fractional period error is small 
outside of an exponentially narrow region in h close to 
AM. 

The potentials reproduce general relativistic precession 



with varying degrees of accuracy and simplicity. Namely, 
these potentials produce accurate relativistic precession: 
(Potential A) in the far field (equation [T3l) ; (Potential 
B) in the far field and with the logarithmic divergence as 
h — )• AM (equation [22]); and (Potential C) in the far field 
and with the correct magnitude of logarithmic divergence 
as h^ AM (equation 1211). 

The potentials described do not include the effects 
of spin, or gravitational radiation, which can be astro- 
physically important for orbits passing close to the hole. 
Neither of these effects can be described by a pseudo- 
Newtonian potential without the presence of undesirable 
inseparable terms including both r and v, and so were 
not considered in this work. For objects whose mass 
ratio with the black hole is sufficiently close to the test 
particle limit, the effects of gravitational radiation can be 
included by subtracting energy and angular momentum 
at per iapsis, for example using the results of I G air et al.l 
(|200l . 

Close to the black hole these potentials should be in- 
terpreted with care since although they diverge at the 
correct angular momentum, the r at which this occurs 
does not correspond to the Schwarzschild radial coordi- 
nate. 
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